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PREFACE 


This  report  documents  the  results  of  USAFETAC  Project  90424,  "Density  Altitude  Maps."  Analysis  were  Major 
Waller  F.  Miller  (DNO)  and  MSgl  Gary  Tryon  (ECO).  The  original  request  (from  OL-A,  Detachment  21,  5  WS, 
Fort  Stewart,  GA)  was  for  maps  of  Iran  and  Iraq  with  coniour.s  of  mcim  monthly  density  altitude  at  the  surface  near 
the  limes  of  maximum  and  minimum  temperature.  Surface  values  of  lempcrdlure.  vapor  pressure,  and  pressure  used 
to  calculate  DA  were  adjusted  at  each  grid  point  using  gridded  terrain  elevation.  The  terrain  used  in  this  analysis 
was  obtained  from  the  Defense  Mapping  Agency  (DMA)  lOO-meter  Digital  Terrain  Database  (DTED).  The 
customer  asked  for  and  received  a  map  scale  of  1:3,(XX),(X)0  in  the  initial  report,  but  a  much  smaller  scale  is  used  in 
this,  the  finhl  publication. 
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1.1  Product  Description.  As  rcqucsicii  by  ihc  cusiomcr,  two  scls  of  planning  maps  arc  provided.  TIic  firsl  set 
provides  contoured  mean  surface  DA  ai  the  time  of  minintimt  temperature,  by  month.  The  second  set  gives  gridded 
monthly  mean  surface  DA  at  the  lime  of  niaximiun  temperature.  The  customer  (Del  21 , 5WS)  asked  for  and  received 
this  data  on  I:3,(XX),000-scaIc  maps  (21.9  inches  x  27.7  inches),  with  contours  at  I,(XX)-fooi  inicr\'als.  These  maps 
are  reproduced  in  the  appendix  of  this  report,  but  to  a  much  smaller  scale.  North-south  map  boundaries  arc  40  and 
25°  N;  east-west  boundaries,  63  and  42°  E. 

1.2  Density  Altitude  Defined.  Density  altitude  (DA)  is  pressure  altitude  (PA)  corrected  for  nonstandard 
temperature  and  dew-point  variations.  PA  is  the  altitude  indicated  by  an  altimeter  when  it  is  set  to  the  standard 
datum  plane  of  29.92  inches  of  mercury.  When  conditions  arc  standard,  PA  and  DA  arc  the  same.  If  temperature  is 
above  standard,  DA  is  higher  than  PA.  If  temperature  is  below  standard,  DA  is  lower  than  PA.  DA  can  also  be 
expressed  as  the  altitude  in  the  ICAO  standard  atmosphere  at  which  a  given  dcasity  (Kciirs.  The  maps  prtxhiccd  in 
this  study  show  the  altitude  in  the  ICAO  standard  atmosphere  at  which  surface  density  occurs.  DA  am  be  calculated 
from  (List,  1984): 


(1) 


where 

DA  =  DA  in  feet 

Psvt  “  station  pressure  in  inches  of  mercury' 

7*^  =  virtual  temperature  in  Rankine  (459-4 +°F) 

Virtual  temperature  can  be  obtained  from  (Dufficld  and  Nasuom,  1983): 

T.  =-0.288+ 1.8(7  +  273.15)  [' ' 

where 

T  =  temperature  in  degrees  C 
r  =  mixing  ratio  (kg/kg) 

Mixing  ratio  is  calculated  by  (Dufficld  and  Nasirom,  1983): 

^_0.622e 

Pstn  "C 


(2) 


(3) 


where 


station  pressure  in  millibars 
vapor  pressure  in  millibars 
mixing  ratio  (kg/kg) 


1 


iL 


A  modified  Tcicn’s  formula  (Murray,  1967)  is  used  to  calculate  vapor  pressure: 


c  =  6.11  {mb)  *  exp  (4) 

where 

r^=  dew  point  icmpcraturc  in 

e  =  vapor  pressure  in  millibars 
a  =  17.269whcnr^>0''Cand2l.874whcnr^<0®C 
b  =  2373  when  7^  >  O'"  C  and  265.5  when  7'^  0®  C 


In  this  study,  the  surface  values  of  temperature,  vapor  pressure,  and  pressure  used  to  calculate  DA  were  adjusted  at 
each  grid  point  asing  gridded  terrain  elevation. 


1.3  Limitations.  The  original  maps  prepared  with  terrain  data  at  G-NM  resolution  were  hard  to  read;  as  a  result, 
data  at  30-NM  resolution  was  interpolated  and  smoothed  to  a  6-NM  grid.  Elevation  errors  of  100  to  700  meters  and 
DA  errors  of  2,000  feet  arc  found  in  areas  of  steep  terrain;  for  this  reason,  the  maps  should  be  ased  for  planning 
purposes  only.  Users  should  also  note  that  DA  in  the  vicinity  of  mountain  peaks  is  significantly  higher  than  the 
surface  DAs  shown  on  the  maps. 

1.4  Data.  Enough  data  was  available  for  38  weather  repoaing  stations  either  in  the  area  of  interest  or  within  1 
degree  of  its  boundaries.  The  terrain  used  for  lliis  analysis  was  taken  from  the  lOO-meter  I>cfcnse  Mapping  Agency 
(DMA)  Digital  Terrain  Database(DTED).  Processing  of  the  DTED  database  is  described  in  Section  3. 

1  Terrain  Effects.  The  variables  from  which  DA  is  calculatcd-tcmpcralurc,  moisture  (mixing  rauo  that  can  be 
determined  from  dew-point  temperature),  and  pressure-are  infiucnced  by  terrain  height,  solar  radiation,  vegetation 
type,  precipitation,  and  geography.  Unfortunately,  there  has  not  been  enough  research  to  take  all  these  factors  into 
consideration;  this  study  only  considers  the  effects  of  terrain,  as  explained  in  the  next  paragraph.  Other  variables 
can  only  be  accounted  for  by  the  representativeness  of  the  observation;  from  available  w^eather  stations. 

1.6  Procedures.  Monthly  means  of  temperature,  vapor  pressure,  and  aliimcicr  setting  were  calculated  for  each 
station  using  values  at  the  lime  of  each  ifaily  maximum  and  minimum  temperature.  These  monthly  mean  values 
were  then  adjusted  to  sea  level  (zero  meters  elevation)  using  the  relationships  described  in  Section  3.  The  \*;triables 
adjusted  to  sea  level  (temperature,  vapor  pressure,  and  altimeter  selling)  were  then  interpolated  to  a  tenth  of  a  degree 
tnuiudc/Iongiiudc  (6-NM)  grid.  This  horizontal  interpolation  was  {Krformal  using  a  cubic  spline  technique  that  fils 
a  third-order  ctiuaiion  in  the  x  (E-W)  and  y  (N-S)  direction.  It  was  now  possible  to  obtain  values  for  c;K:h  grid  point 
at  sea  level.  The  gridded  sea-level  meteorological  data  and  the  gridded  terrain  heights  were  then  used  to  adjust 
temperature  and  vapor  pressutc  to  tlie  terrain  surface  using  suindard  meteorological  relationships.  Altimeter  .setting 
was  adjusted  to  station  (surface)  pressure  using  the  gridded  elevation.  Once  these  variables  were  obiaincd  at  the 
surface.  Equation  I  could  be  used  to  calculate  DA.  These  steps  arc  described  in  greater  detail  in  Section  3. 
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1.7  QualKy  Control.  The  following  assumpiions  were  necessary  to  complete  ihLs  study: 

'That  DA  at  the  time  of  the  daily  maximum  temperature  is  close  to  the  maximum  DA  for  the  day. 

'That  minimum  DA  has  to  be  related  to  DA  at  the  minimum  temperature. 

'That  the  results  would  net  be  seriously  affected  by  the  lack  of  available  data  for  Iran  and  Iraq-that  is, 
the  number  of  limited  duty  stations  and  3-houriy  observations,  and  a  period  of  record  less  than  10  years. 

'That  the  cubic  spline  interpolation  .spreads  the  data  accurately. 

The  error  introduced  by  each  of  these  assumptions  will  be  discu.sscd  in  Section  4,  where  support  for  the  complex 
procedure  used  in  producing  our  DA  maps,  as  oppo.scd  to  interpolating  DA  directly  to  a  grid,  is  also  supplied. 


1.8  Maps.  The  D.A  maps  arc  described  in  Section  5  and  provided  in  the  appendix.  DA  isopleths  arc  labeled  at 
I.OOOTool  intervals.  Because  of  the  rugged  terrain  in  Iran,  .smoothed  terrain  was  used;  this  impro%'es  readability  at 
the  expense  of  some  accuracy.  To  take  full  advantage  of  the  high-resolution  terrain  data  available  from  DMAAC, 
much  larger  maps  would  be  required. 
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2.  INPUT  DATA 


2.1  Surface  Observations.  Alihougli  USAFETAC  mainiains  surface  observations  from  1973  lo  1989,  some 
stations  reported  only  at  3-hourIy  intervals  and  others  have  failed  to  report  for  several  years.  Rodney  (1986) 
concluded  that  10  years  of  observational  data  was  enough  for  a  representative  sample  of  climate  in  a  region.  Only 
1 1  of  the  38  stations  used  here,  however,  provided  a  full  10*year  period  of  record  (FOR).  Most  stations  had  a  5-ycar 
FOR,  but  stations  with  FORs  as  small  as  4  years  were  included  if  they  w'cre  in  a  data-spaise  area.  Table  1  lists  the 
weather  stations  used;  their  locadons  are  shown  in  Figure  1 .  Some  stations,  although  olT  the  maps,  were  included  as 
inputs  to  the  interpolation  scheme  to  provide  mere  accurate  values  at  the  boundaries.  As  can  be  seen  in  Rgure  I, 
most  of  tile  stations  are  in  the  western  half  of  tlic  region,  while  the  eastern  half  has  very  few.  We  evaluated  35 
stations  in  the  eastern  portion  and  could  find  none  with  enough  observations  to  be  useful. 


TABLE  I.  List  of  Stations  Used  in  the  Interpolation  Scheme. 


Bk>ck 

Station 

Nam*  Country 

Lat 

Lon 

El*v  (meters) 

FOR 

170965 

Erzunixn 

TU 

3957N 

41  lOE 

1756 

77.86 

172020 

Ela/ig 

TU 

3836’N 

39  17-E 

882 

77.86 

172800 

Do'atbakir 

TU 

3703’N 

40irE 

677 

77.86 

379070 

Rzuli 

RS 

3936-N 

4709*E 

430 

73.77 

379850 

Lenkonn 

RS 

3844-N 

4S50*E 

.11 

73.77 

385070 

Kr25nov0<I^ 

RA 

4002’N 

5259T 

S9 

73.77 

388SOO 

Ashabad 

RA 

375S-N 

5820E 

228 

73.77 

403570 

Am 

SD 

3054'N 

4lt^-E 

552 

77.86 

403620 

Kstha 

SD 

293S'N 

43  ^’E 

447 

77.86 

403730 

HzfrATBatln 

SD 

282D-N 

4607’E 

355 

77.86 

403940 

Hril 

SD 

2726*N 

414VB 

1013 

77.86 

404050 

Gatssim 

SD 

26I$'N 

4346*E 

6^ 

77.86 

404160 

Dhahrsn  lAP 

SD 

26 16'N 

5009E 

17 

77.86 

401380 

Riyadh 

SD 

244rN 

4643-E 

612 

77.86 

405S20 

Kuwait  !AP 

KW 

2913’N 

4759'E 

55 

77.86 

406110 

Sahluddin 

IQ 

3637'N 

441TE 

loss 

77.80 

406210 

Kirkuk 

IQ 

352S’N 

4424*E 

331 

73.80 

406230 

Suhimanha 

IQ 

3533N 

452TE 

853 

7340 

406290 

Arm 

IQ 

342S-N 

41 57'E 

150 

7340 

406340 

Hzdiiha 

IQ 

3404'N 

4222E 

140 

7340 

406370 

Kanaqci 

IQ 

34  IS’N 

4526‘E 

202 

7340 

406120 

Rulbah 

IQ 

330rN 

40 17-E 

615 

7340 

4D650P 

B2|;|K!2d 

IQ 

33  i4-N 

4414*E 

34 

7340 

406560 

Ksxfulaa 

IQ 

3237N 

44  ore 

7740 

4!^i50 

Kmaldia* 

IQ 

32  iO’N 

4603*E 

15 

7340 

406700 

IQ 

31 59  N 

44  I9*£ 

32 

7340 

406720 

Diw^ntya 

IQ 

31 59’N 

4.559’E 

20 

7340 

406740 

Semrars 

IQ 

3!  IS  N 

45 16’E 

6 

7340 

406760 

Naririya 

IQ 

31 05’N 

46  i4*E 

3 

7340 

4mm 

Amzndi 

IQ 

315rN 

47  lO’E 

9 

7340 

406S90 

Rasnh/Mzgzl 

IQ 

3034-N 

4747*E 

7340 

407360 

Babubsr 

IR 

3643’N 

523rE 

.21 

7340 

407450 

Mashhad 

IR 

36  i6-N 

593S'E 

9S0 

7640 

407650 

TchnnfMchrrbad 

IR 

3.54rN 

51  are 

1191 

7340 

407690 

Ar^ 

IR 

3406*N 

494Z-E 

1729 

7640 

408030 

Esfahan 

IR 

3237-N 

5140'E 

1590 

7640 

408310 

Ab^ianfAr 

IR 

3022‘N 

48 15’E 

11 

7340 

41 1950 

Shaijah 

ER 

2520‘N 

553rE 

33 

7746 
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2:2  Variables  Used.  For  this  sludy,  icmpcralure,  dcw-poinl  icmpcraliire,  and  allimelcr  selling  were  used  lo 
calculalc  DA  (see  Equalion  1  and  1.6).  U  allimelcr  selling  was  nol  available,  sea-level  pressure  was  convened  lo 
allimelcr  selling  as  described  In  Section  3.  Allimelcr  selling  was  chosen  for  three  reasons: 

•Station  pressure  is  nol  reponed  in  surface  observations. 

•Obtaining  station  pressure  from  altimeter  selling  o.  .y  requires  one  additional  variabic-cicvalion. 

•The  technique  used  in  this  sludy  requires  that  pressure  be  adjusted  lo  sea  level  (zero  meters)  for 

interpolation. 

Vapor  pressure  was  obtained  from  Equation  4  using  ihc  dcw-poinl  icmperalure.  Monthly  means  of  temperature, 
vapor  pressure,  and  allimelcr  selling  were  calculated  using  the  values  at  the  limes  of  the  maximum  and  minimum 
icmperalure.  Some  variables  are  in  metric  units  and  some  in  English  units,  depending  on  which  calculation  was 
being  done.  The  equations  in  this  repori  are  in  the  form  given  in  the  reference  from  which  ihcy  were  taken; 
therefore,  all  equations  specify  units  required. 

2.3  Terrain  Data.  Although  USAFETAC  has  arranged  lo  acquire  worldwide  DTED  terrain  data  at  lOO-mclcr 
resolution  from  the  Defense  Mapping  Agency  (DMA),  only  a  small  subscl  of  the  DTED  database  was  available  lo 
USAFETAC  analysts  for  use  in  this  sludy.  Fortunately,  the  daia  that  was  available  included  Iran  and  Iraq.  To 
conserve  storage,  only  enough  data  arc  maintained  at  USAFETAC  to  recreate  the  actual  data  with  less  than  1  percent 
of  interpolated  values  having  errors  greater  than  3(X)  feet.  Terrain  heights  with  6-NM  resolution  were  obtained  from 
this  dataset  for  the  area  of  interest;  as  shown  in  Figure  2,  using  terrain  data  to  this  resolution  would  result  in  an 
unreadable,  unusable  product 
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Figure  2.  FJevations  (meters)  for  Iran  and  Iraq  Using  the  DTFI)  Database  with  6-NM  Resolution. 


3.  PROCEDURES  FOR  DA  CALCULATION 


3.1  Gross  Error  Checks,  Surface  weallier  observations  for  each  station  were  quality  controlled  for  gross  errors. 
Tlie  data  had  to  pass  the  following  tests: 

•Temperature  must  be  between  20®  and  138®  F. 

•Dew-point  lempcralure  must  be  less  Ilian  or  equal  to  the  temperature. 

•Dew-point  temperature  must  be  between  -30®  and  93®  F. 

•If  altimeter  setting  is  used,  it  must  be  between  28  and  32  inches  of  mercury. 

•If  sea  level  piessure  is  used,  it  must  be  between  948  and  1050  mb. 

•The  difference  between  the  maximum  and  minimum  temperature  must  be  less  than  65®  F. 

•The  difference  between  the  altimeter  setting  for  the  minimum  and  maximum  temperature  is  limited  to 
half  an  inch  of  mercury;  the  sea-level  pressure  difference  is  limited  to  20  mb. 

If  all  these  conditions  were  met,  the  surface  weather  observation  was  used  in  selecting  the  maximum  and  minimum 
temperature  as  described  in  3.2.  Extreme  highs  for  temperature  and  dew-point  temperature  used  in 
quality-controlling  the  data  were  obtained  from  MIL-STD-210C.  The  lower  values  were  thought  to  be  appropriate 
extremes  for  the  Mideast. 

3.2  Maximum  and  Minimum  Values.  The  maximum  and  minimum  temperatures  for  each  day  were  selected, 
along  with  the  dew-point  temperature  and  altimeter  setting  (sea  level  pressure)  that  corresponded  to  the  time  of  the 
maximum  and  minimum  temperature.  Because  the  data  from  some  stations  was  limited,  only  minimum  and 
maximum  temperatures  occurring  within  3  hours  of  0300  and  1200Z  (respectively)  were  selected. 

3.3  Vapor  Pressure.  Vapor  pressure  was  calculated  using  the  dew-point  temperature  at  the  lime  of  the 
maximum  and  minimum  temperature  for  each  day  (Equation  4).  The  monthly  mean  vap(>r  pressure  was  then 
calculated  for  these  two  times.  Another  option  would  have  been  to  calculate  monthly  mean  dew-point  temperatures 
at  the  lime  of  maximum  and  minimum  temperatures  and  use  those  values  in  Equation  4  to  calculate  the  monthly 
mean  vapor  pressure.  Tests  showed  that  there  was  a  difference  as  large  as  0.5  mb  when  the  monthly  mean  dew-point 
temperature  was  used;  we  therefore  calculated  monthly  mean  from  daily  vapor  pressure  at  the  time  of  the  maximum 
and  minimum  temperature. 

3.4  Pressure  Conversion.  Altimeter  setting  was  chosen  for  interpolation  to  the  horixonial  grid  because  only 
terrain  heights  were  needed  to  obtain  surface  pressure.  Several  stations  used  in  this  study,  however,  reported  only 
sea  level  pressure;  for  them,  wc  converted  sca-K pressure  to  altimeter  setting  by  reversing  the  sea-level 
computation  process  given  in  llie  Smithsonian  Meteorological  Tables  (List,  1984).  Mean  surface  temperature  was 
calculated  from  daily  minimum  and  maximum  temperature.  To  calculate  mean  layer  temperature,  surface 
temperature  was  increased  by  5®  C  for  every  1,000  meters  the  station  was  above  sea  level  (List,  1984).  Since  the 
mean  layer  temperature  is  valid  hall  the  distance  between  the  surface  and  sett  level,  the  lapse  rale  term  can  be  added 
to  the  mean  surface  temperature  if  it  is  reduced  by  half;  llial  is,  (0.5(5®  C)/l  ,(XK)  meters  or  l/4(K)  ®C  m'*).  The  mean 
temperature  for  the  layer  is  used  to  solve  the  hypsomcu^ic  equation  for  station  pressure  (Wallace  and  Hobbs,  1977): 


Pstn  =Psl  exp 


(-g  ft) 


(5) 


8 


where 


Pstn  ~  pressure  in  mb 
p^j  =  sea-level  pressure  in  mb 

g  =  acceleration  due  to  gravity  in  ms'^ 
h  =  sUttion  elevation  in  meters 
R  =  gas  constant  for  dry  air,  287  mV^K'* 

=  mean  layer  temperature  in  kelvin  (see  section  3.4),  given  by 

'L  =  0*5  (T^  -h  7' . J  +  hm) 

^max  "  ni^ixinium  temperature  for  the  day  in  kelvin 
^min  “  *if^inimum  temperature  for  the  day  in  kelvin 

Station  pressure  is  than  converted  back  to  altimeter  setting  using  Equation  6  (Duffield  and  Nastrom,  1983): 

ALSTG  =  (1  +  8.42  •  k  ps,n  ^  Ps,n  ^§§5  (6) 


where 


ALSTG  =  altimeter  setting  in  inches  of  mercury 
Pstn  “  station  pressure  in  mb 


3.5  Monthly  Mean.  The  monthly  mean  of  temperature,  vapor  pressure,  and  altimeter  setting  is  calculated  using 
the  tcmjKiralurc,  vapor  pressure,  and  altimeter  setting  at  the  time  of  maximum  temperature  for  each  day.  This  time 
must  be  within  3  hours  of  120()Z,  but  it  docs  not  have  to  be  the  .same  time  every  day.  This  step  is  repeated  for  values 
at  the  time  of  the  daily  minimum  temperature,  which  is  within  3  hours  of  03(K)Z.  The  standard  lapse  rate  of  6.5®  C 
per  l,(X)0  meters  is  used  to  adjust  the  monthly  mean  temperature  to  .sea  level  (List,  1984): 

Tsi^Tsfc^O.mSh  (7) 


where 

k  =  station  elevation  in  meters 
Tji  5=  temperature  at  sea  level  in  ®C 
7'^^  =r  temperature  at  the  surface  in  ®C 
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Hann’s  empirical  formula  a?  described  by  List,  1984,  is  used  lo  adjust  the  monthly  mean  vapor  pressure  to  sea  level; 
altimeter  setting  is  already  adjusted  lo  a  constant  level. 


Csi  =  Csfc  exp 


(8) 


where 

=  vapor  pressure  at  sea  level  in  mb 
Vc  ~  vapor  pressure  at  the  surface  in  mb 


3.6  Interpolation.  Sea-level  temperature,  sea-level  vapor  pressure,  and  altimeter  setting  at  irregularly  spaced 
stations  are  interpolated  to  a  regularly  spaced,  horizontal  (scit-Ievel)  grid  using  a  cubic  spline  technique.  The  grid 
consists  of  points  every  tenth  of  a  degree  in  latitude  and  longitude.  Examples  of  interpolated  data  for  July  arc 
provided  in  Figures  3  through  5.  A  cubic  spline  produces  a  surface  similar  lo  "one  that  would  be  formed  if  a  stiff, 
thin  metal  plate  were  forced  through  or  near  the  given  data  points"  (SAS,  1990).  Mathematically,  a  cubic  spline  fils 
a  third-order  polynomial  lo  the  data  in  both  horizontal  directions.  Values  for  each  grid  point  can  be  determined  from 
the  third-order  equations.  Thus,  data  for  38  irregularly  spaced  points  arc  spread  lo  each  point  on  the  tenth  of  a 
degree  grid  at  sea  level.  Extension  of  data  in  the  vertical  and  compulation  of  DA  is  a  subsequent  processing  step 
and  is  described  in  3.8,  below. 

3.7  Smoothed  Terrain,  As  this  report  was  l)cing  publi.shcd,  USAFETAC  was  still  in  the  process  of  acquiring 
DTED  data  from  DMA.  To  mainuiin  a  DTED  database  for  the  entire  world,  USAFETAC  would  have  had  to  find 
room  to  store  6,300  magnetic  tapes.  A  mclhtxi  was  developed,  however,  to  store  just  enough  daut  to  rcprwlucc  the 
original  terrain  with  less  than  1  percent  of  die  errors  greater  than  3(K)  feet.  From  this  reduced  data  set,  we  obtained 
6-NM  resolution  terrain  for  the  area  of  interest.  Initial  studies  using  this  data  showed  that  there  was  too  tnuch  detail 
"SCO  Figure  2.  To  provide  smoother  fields,  30-nm  terrain  data  was  intcr|)olalcd  using  the  .same  cubic  spline 
technique  as  for  the  meteorological  variables.  During  interpolation,  the  data  was  smoothed;  the  smoothed  terrain 
field  is  shown  in  Figure  6.  A  comparison  of  Figures  2  and  6  shows  the  detail  lost;  the  smoothed  map,  however,  is 
much  more  readable. 

3.8  DA  Calculations.  After  interpolation,  we  had  monthly  mean  .sca-lcvel  temperature,  monthly  mean  sca-lcvcl 
vapor  pressure,  monthly  mean  altimeter  setting,  and  elevation  for  each  grid  point.  There  are  six  grids  for  each 
month;  one  set  of  three  (sca-levcl  temperature,  sca-lcvel  vapor  pressure,  and  altimeter  setting)  for  the  time  of 
maximum  temperature  and  another  .set  of  three  for  the  lime  of  minimum  temperature.  The  procedure  for  calculating 
DA  from  gridded  fields  was  repeated  for  citch  month  and  for  each  set  of  grids.  Equations  7  and  8  arc  solved  for 
.surface  temperature  and  vapor  pre.ssure  using  the  smoothed,  gridded  terrain  heights.  Station  prc.ssurc  is  obtained 
from  altimeter  setting  using  the  rclalion.ship  from  Dufficid  and  Nasirom,  1983  (Equation  9);  finally,  DA  was 
calculated  with  Equation  1. 


!),,„=  {ALSTG  -4.3()6-  l(r* 


(9) 


where 

ALKIG  =  altimeter  selling  in  inches  of  mercury 
Pstn  ”  station  pressure  in  inches  of  mercury 
h  =  gridded  terrain  elevation 
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Figure  4.  Interpolated  Mean  Monthly  Sca-Ixvcl  Vapor  Pressure  (mb)  at  the  Time  of 
Maximum  Temperature,  July. 
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Figure  5.  Interpolated  Mean  Monthly  Altimeter  Settings  (in  Hg)  at  the  Time  of  Maximum  Temperature,  July. 
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Figure  6.  Elevations  (meters)  from  Smoothed  Terrain  Data  for  Iran  and  Iraq. 
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Figure?,  Interpolated  Mean  Monthly  Maximum  Sea^Level  Temperature  (^C),  July;  12  Stations  Removed, 


Figure  8,  Mean  Monthly  Surfac<*  Density  AKidude  (feet)  at  the  Time  of  Maximum 
T€mpcra(ure,July;  Terrain  Effect:;  inelUucd^ 


Figure  P.  Mean  Monthly  Surface  Density  Altitude  (feet)  at  the  Time  of  ^faximum 
Temperature,  July;  Terrain  Effects  included,  but  with  12  Stations  Removed. 


Figure  10.  Interpolated  Mean  Monthly  Surface  Density  Altitude  (feet)  at  the  Time  ai' 
Maximum  Temperature ,  July;  Terrain  Effects  Not  Included. 


4.  QUALITY  CONTROL 

4.1  Accuracies  of  Maximum  and  Minimum  Temperatures  were  checked  against  published  valuas  in 
USAFETAC/DS-8yA)35,  Station  Climatic  Summaries  Asia.  More  ilum  half  the  suuions  used  could  be  verified  by 
this  method.  Errors  were  on  the  order  of  2  to  3®  F.  Even  though  only  4  years  of  data  were  available  for  some 
stations^  results  from  these  small  samples  seemed  to  be  reprasentative  of  the  climatological  mean. 

4.2  DA  at  Maximum  Temperature.  Our  gridding  mciluxl  used  data  valid  at  die  time  of  die  daily  maximum  and 
minimum  temperature  to  calculate  DA.  This  method  agrees  widi  die  customer’s  request  for  DA  at  the  lime  of  the 
maximum  and  minimum  jcmpcralurc.  Our  first  quality  control  check  was  to  determine  the  size  of  the  difference 
between  the  monthly  mean  maximum  DA  and  the  monthly  mean  DA  at  the  lime  of  maximum  daily  temperature. 
For  the  six  stations  evaluated,  the  monthly  mean  maximum  DA  was  no  more  Uian  100  feel  higher  than  the  monthly 
means  using  the  maximum  temperature.  An  error  of  similar  size  (opposite  sign)  is  associated  with  the  monthly 
mean  minimum  DA  and  the  monthly  mean  D.A  at  the  lime  of  the  minimum  temperature. 

4.3  Mean  DA  From  Mean  Values.  The  next  concern  was  if  the  DA  calculated  from  die  monthly  mean  values  of 
temperature,  vapor  pressure,  and  altimeter  setting  was  close  to  the  monthly  mean  maximum  and  minimum  DA.  The 
error  belwech  the  monthly  mean  DA  and  die  calculated  mean  was  100  feel.  Our  results  showed  excellent  agreement 
with  DA  calculated  from  the  values  at  the  time  of  maximum  or  minimum  temperature  and  the  calculated  monthly 
mean  from  these  values.  The  errors  were  lass  dian  20  feet. 

4.4  Thrce^Kourly  Data.  Next,  we  ran  a  test  using  four  suiiions  to  compare  monthly  mean  DA  at  the  lime  of 
maximum  oj  minimum  temperature  from  hourly  and  3-hourIy  data.  The  difference  between  the  monthly  means 
changed  by  less  than  1(X)  feet  (maximum,  lower;  minimum,  higher). 

4.5  Interpolation  Error.  Because  data  from  several  stations  was  used  to  determine  die  value  at  each  grid  point,  a 
small  error  was  expected.  If  no  data  was  present,  the  size  of  the  error  was  unknown.  Because  of  the  limited  number 
of  stations  in  the  region,  the  enor  at  a  large  number  of  grid  points  was  unknown.  In  an  attempt  at  estimation,  12 
stations  widi  data  were  removed  from  die  interpolation.  Table  2  shows  the  difference  between  interpolated  valuas 
using  all  stations  and  the  inicrpt -ated  value  missing  the  12  stations  (All  -  Void)  for  sea-level  temperature,  .sea-level 
vapor  pressure,  altimeter  setting,  and  DA.  As  might  be  expected,  die  errors  were  larger  at  the  eastern  locations  (the 
last  three  stations).  Figure  7  showeti  die  interpolated  maximum  sea  level  temperature  for  July.  A  comparison  of 
Figures  3  and  7  shows  that  the  basic  pattern  remains,  but  dial  some  of  the  small-scale  features  arc  lost.  These 
differences  are  even  smaller  for  the  monthly  mean  DA  at  the  time  of  maximum  temperature  in  July-sec  Figuras  8 
and  9. 

4.6  Terrain  Error.  The  largest  source  of  error  was  the  elevation  assigned  to  each  grid  point.  Because  of  terrain 
resolution  and  smoothing,  there  were  diffcrcncas  of  up  to  7{)()  meters  in  the  comparison.  Elevation  errors  this  large 
can  result  in  DA  errors  on  the  order  of  a  2,000  feel.  When  the  gridded  elevation  was  close  to  station  elevation,  the 
DA  was  within  several  hundred  feet.  This  error,  however,  wiis  a  compromise  between  accurate  data  and  readable 
data. 

4.7  Spreading  DA.  One  might  ask,  "Why  go  llirough  the  complex  procedure  of  taking  elevation  into  account 
when  creating  gridded  fields  of  DA?"  A  comparison  of  Figuras  10  (interpolated  DA)  and  8  (terrain  effects  included) 
shows  the  large  errors  that  rasult  from  ignoring  terrain,  ‘fliase  errors  were  quantified  using  die  daui  void  last  rasulLs. 
Table  2  shows  the  difference  between  the  mean  monthly  DA  at  the  lime  of  die  maximum  icmp'ralurc  and  the 
interpolated  DA  (Spread  DA).  The  mean  monthly  DA  calculated  witli  terrain  (but  with  12  missing  Ions)  is  also 
compiircd  to  the  station  value  O'^ble  2,  Terrain  DA).  The  groiiast  errors  occur  in  area  with  the  highasi  terrain. 

Data  was  taken  from  DA  at  the  lime  of  maximum  temperature  for  July.  Spread  DA  is  die  diffc'  :;t;cc  between  station 
DA  and  interpolated  DA  widiout  taking  terrain  into  acouni,  and  with  12  .stations  removed.  Terrain  DA  is  the 
difference  between  station  DA  and  DA  calculated  with  terrain,  but  v/iih  12  stations  removed. 


15 


TABLE  2.  DifTerences  Between  Selected  Interpolated  Variables:  Sprea.^  DA  and  Terrain  DA  (feet)* 


BLKSTN 

o 

O 

e„(mb) 

ALSTG  (In  hg) 

Spread  DA  (feet) 

Terrain  DA  (feet) 

385070 

3.6 

-3.9 

-0.05 

-331 

791 

403620 

0.5 

0.4 

-0.01, 

345 

212 

403730 

-0.1 

0.1 

0.03 

456 

130 

406110 

0.7 

0.6 

-0.10 

-536 

292 

406210 

0.5 

0.8 

0.03 

1,089 

186 

406700 

-0.3 

-2.6 

0.00 

96 

96 

406890 

-2.4 

8.1 

-0.01 

163 

20 

407540 

12.0 

-9.8 

0.10 

-4,615 

1,(18 

407690 

7.8 

-6.2 

0.26 

-5,630 

153 

408000 

7.2 

-8.7 

0.19 

-5,117 

715 

4.8  Total  Error.  Based  on  ihc  error  analysis  above  and  ignoring  ihc  smoothing  of  terrain  height,  the  DAs 
provided  arc  accurate  to  within  1  ,(XX)  feet.  The  errors  along  the  western  half  of  the  study  area  arc  less  than  5(X)  feet 
because  of  the  higher  data  density.  Smoothing  of  terrain  can  result  in  another  1,000-foot  error.  The  DA  at  tops 
of  large  mountains  is  not  considered  in  this  analysis. 


16 


5.  MAPS 


5.1  Description.  Twenly-four  maps  of  mean  DA  for  Iran  and  Iraq  are  provided  in  the  appendix.  The  map  area, 
specified  by  the  customer,  is  bounded  by  25  and  40®  N  latitude  and  43  and  62®  E  longitude.  There  is  a  map  for  each 
month  at  the  time  of  maximum  and  minimum  temperature.  Contours  are  drawn  every  1,000  feet,  with  zero  feet  as 
the  base.  Table  3  provides  DA  values  for  Baghdad,  Iraq;  Tehran,  Iran;  and  Dhahran,  Saudia  Arabia,  for  use  as 
reference  points.  Geopolitical  boundaries  of  the  area  (bold  lines)  arc  also  plotted.  Our  software  docs  not  close  off 
southern  Iraq  at  the  Persian  Gulf  coast.  Latitude  and  longitude  lines  arc  provided  at  5-dcgrcc  increments  to  help 
users  locate  points. 

5.2  Map  Data  Summary.  Minimum  DAs  are  found  in  Iraq,  Saudi  Arabia,  and  Russia,  but  DA  increases  rapidly 
along  the  border  of  Iran,  reaching  its  largest  values  on  the  central  plateau.  During  the  winter,  the  DA  at  the 
minimum  temperature  along  the  Persian  Gulf  is  near  -1,000  feet,  increasing  to  several  thousand  feet  during  the 
summer.  DA  at  the  maximum  temperature  is  near  zero  during  die  winter,  increasing  to  almost  4,000  feet  during  the 
summer.  In  the  high  plateau  of  Iran,  DA  at  the  minimum  temperature  is  near  5,000  feet  during  the  winter  and  7,000 
feet  during  the  summer.  DA  at  the  maximum  temperature  is  more  than  8,000  feet  during  the  winter,  increasing  to 
more  than  1 1,000  feet  during  the  summer. 


TABLE  3.  DA  (feet)  for  Key  Reference  Point*!. 


Month 

Baghdad 

Minimum  Maximum 

Tehran 

Minimum  Maximum 

Dhahran 

Minimum  Maximum 

January 

-1876 

-363 

2706 

5741 

-682 

560 

February 

-1451 

201 

3173 

6170 

-471 

808 

March 

-905 

776 

3885 

6926 

80 

2004 

April 

-170 

1570 

4763 

7784 

798 

2121 

May 

-64 

2421 

5370 

8387 

1489 

2873 

June 

1180 

3114 

6135 

9196 

2064 

3419 

July 

1548 

3526 

6554 

9573 

2350 

3703 

August 

1360 

3374 

6343 

9302 

2241 

3546 

September 

789 

2910 

5775 

8768 

1714 

3153 

October 

-137 

1995 

4790 

7879 

1039 

2417 

November 

-1194 

765 

3741 

6791 

280 

1541 

December 

-1602 

-68 

3190 

6140 

-466 

793 
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6.  SUMMARY 


6.1  Terrain  Effects.  Mean  DA  for  Iran  and  Iraq  is  provided  for  each  monih  al  ihe  lime  of  maximum  and 
minimum  icmperaiurc.  Instead  of  inicrpolaling  mean  DA  for  all  38  available  stations,  terrain  elevation  and  standard 
meteorological  relationships  were  used  to  provide  DA  on  a  grid.  This  technique  required  converting  temperature 
and  vapor  pressure  to  a  common  surface  (sea  level;  before  interpolating  them  to  a  laiiludc-Iongitudc  grid  with  a 
resolution  of  O.I  degree.  These  interpolated  values  were  than  adjusted  back  to  the  gridded  elevation.  From  these, 
adjusted  values  of  temperature,  vapor  prc.ssure,  and  pressure,  DA  were  calculated.  Much  more  detail,  along  with 
lighter  gradients  of  DA.  were  produced  by  this  method  than  by  merely  interpolating  mean  monthly  DA. 

6.2  Quality.  An  error  analysis  .showed  that  the  DAs  provided  have  errors  of  less  than  500  feet  in  the  western  half 
of  the  region,  increasing  to  1,000  feel  in  the  east.  Additional  error  is  added  by  using  smooth  terrain  instead  of  the 
full  6-NM  re.solulion;  however,  the  smooth  terrain  was  justified  in  that  its  use  resulted  in  a  readable  producL 


APPENDIX 


The  24  maps  in  the  appendix  provide  contoured  values  of  density  altitude  (DA)  in  Iraq  and  Iran  at  minimum  and 
maximum  temperature  for  each  month  of  the  year.  Values  of  DA  were  calculatai  as  described  in  the  text  of  this 
reporL  Contours  were  drawn  with  SAS  .softw-sirc,  a  copyright-iKolcclcd  product  of  the  SAS  Institute.  Caiy,  NC 
Contours  arc  in  l.OOOTooi  intcr\*als.  Since  SAS  .software  was  incapable  of  labeling  contours,  the  labels  were 
applied  by  hand. 
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Figure  A-1.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Minimum  Temperature,  January. 


Figure  A*3.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Minimum  Temperature,  February. 


Figure  A-4.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Maximum  Temperature,  February. 


Temperature,  March. 


Figure  A-8.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Maximum  Temperature,  April, 
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Figure  A-10.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Maxlmuni  Temperature,  May. 
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Figure  A-12.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  MaxImumTemperature  for  June. 
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Surface  DA  (feet)  for  Iran  and  Iraq  at  Minimum  Temperature,  August. 


Surface  DA  (feet)  for  Iran  and  Iraq  at  Maximum  Temperature,  August. 
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Figure  A-20.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Maximum  Temperature,  October. 
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Figure  A-22.  Mean  Surface  DA  (feet)  for  Iran  and  Iraq  at  Maximum  Temperature,  November. 
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SPECIALIZED  TERMS  AND  ACRINABs 


ACRINAB 

acronym,  inilialistn,  abbrcvhulon 

ALSTG 

altimeter  selling 

C 

Celsius 

DA 

Density  altitude 

DMA 

Defense  Mapping  Agency 

DTED 

Digital  Terrain  Elevation  Database 

e 

vapor  pressure 

ER 

United  Arab  Emirates 

®ji 

sea  level  vapor  pressure 

®5fc 

surface  vapor  pressure 

F 

Fahrenheit 

g 

gravity 

h 

station  elevation 

IQ 

Iraq 

IR 

Iran 

KW 

Kuwait 

mb 

millibar 

P 

pressure 

P,i 

sea-level  pressure 

P*tn 

station  pressure 

POR 

period  of  record 

R 

gas  constant  for  dry  air 

r 

mixing  ratio 

RA 

USSR-ASIA 

RS 

USSR-Europe 

T 

Temperature 

Tj 

Dew-point  temperature 

Mean  layer  temperature 

T 

'max 

Maximum  temperature 

Minimum  temperature 

T„ 

Sea-level  temperature 

T,rc 

Surface  temperature 

Tv 

Virtual  temperature 
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